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We have studied the one-meson radiative tau decays: — > Ur'n:^ {K^)^, computing the 
structure dependent contributions within a Lagrangian approach based on the large-iVc" 
hmit of QCD that ensures the proper low-energy limit dictated by chiral symmetry. Upon 
imposing the short-distance QCD constraints to the form-factors we are able to predict 
the structure dependent radiation without any free parameter and, therefore, the relevant 
observables for the decay — !■ Vt-t^^^: the photon energy spectrum, the invariant mass 
spectrum of the meson-photon system and the integrated decay rate. We also discuss the 
remaining uncertainties in these observables for the t" — > VrK~^ decay. According to our 
results, the present facilities could detect these rare decays for the first time in the near 
future allowing for a stringent test of our predictions. 
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I. INTRODUCTION 

The decays of the r lepton are an ideal benchmark to analyze the hadronization of the spin-one 
QCD currents and to learn about the treatment of the strong interaction involving the intermediate 
meson dynamics in rather clean conditions, since the electroweak part of the decay is under good 
theoretical control and the light-flavoured hadron resonances rule these processes [1-12]. 

Among the various exclusive tau decays, we are particularly interested in the one meson radiative 
tau decays in this article, i.e. r — t- ■k{K)v^. Unlike in many of the multi-pseudoscalar tau decay 
channels, where the excited vector resonances (such as p' , K*') could play a significant role due 
to the kinematics effects, the channels we are going to study here should be less influenced by 
the excited resonances, since the lowest vector resonances such as /3(770) ,i^*(892), could already 
decay into ^{K)^ at their on-shell energy regions. So these radiative decay processes can provide 
us with an excellent environment to investigate the lowest resonance states. On the experimental 
side, these channels have not been detected yet, which makes our analysis more meaningful, since 
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the current work could serve as a helpful tool for the measurements of these channels in the B and 
tau-charm Factories in the near future. 

The decay amplitude for the one meson radiative decay of the r includes an internal 
bremsstrahlung {IB) component, that is given by QED, and thus can be calculated unambiguously 
to any desired order in perturbation theory. In addition, one has the structure dependent (SD) 
part, dominated by the effects of the strong interaction. Lorentz symmetry determines that there 
are two independent structures, the so-called vector and axial-vector form factors that encode our 
lack of knowledge of the precise mechanism responsible for hadronization. 

One has then to rely on parametrizations of these form factors. Some of the earliest attempts can 
be found in Refs. [1]. The so-called Kiihn-Santamaria model [G] became a very popular approach. 
For a given meson mode it proceeds as follows: it normalizes the form factor using the lowest order 
contribution stemming from the low-energy effective field theory of QCD in the light-flavoured 
sector: Chiral Perturbation Theory, X-P^ [13-16]. Then, the form factor is constructed in terms of 
Breit-Wigner functions weighted by some unknown parameters (to be fitted to data) in such a way 
that it vanishes at infinite transfer of momentum in order to obey a Brodsky-Lepage behaviour 
[17, 18]. This procedure was successful to describe the ARGUS data on the r — )■ TTTTirVr decays [19]. 
However, as the data became more precise [20] it was shown that there was room for improvement 
[8, 12]. Among the theoretical reasons for this, one finds that in the low-energy limit this model is 
not consistent with xPT at 0{p^) [8, 21]. Assorted versions of the Kiihn-Santamaria model have 
been used for several two- and three-meson tau decays [3, 22] and also for the radiative one-meson 
decays we consider in this article [4]. In general, these works made following Ref. [•)] suffer from 
additional problems, as discussed in Ref. [23]. In fact, arbitrary parametrizations are of little help 
in order to comprehend the characteristics of the hadronization procedure independently of their 
eventual success in describing experimental data. The best procedure in order to investigate this 
problem should be to build the relevant form factors upon as many QCD features as possible. 

The TAU OLA library [21] is currently the most relevant tool to analyze tau decay data. Al- 
though it only incorporated the Kiihn-Santamaria models at first, it was enriched with parametriza- 
tions by the experimental collaborations ALEPH and CLEO -collected in Ref. [25]- in order to 
explain their data. In the B-factories era, however, it has become evident that this strategy is not 
adequate and TAUOLA has been opened to the introduction of matrix elements obtained within 
other approaches that include more properties stemming from QCD. This makes it an excellent 
tool to exploit the synergies of Monte Carlo methods and theoretical approaches to better under- 
stand the large data samples of high quality obtained at the B-factories (BaBar, Belle) and the 
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r-charm factories, such as BES-111. This is thus an appropriate benchmark where the results of 
our research can be apphed. 

Although one knows the underlying theory, QCD, this kind of study is rather involved since 
there is no fully analytic way of relating the final state mesons that are detected to the quark and 
gluon degrees of freedom of the QCD Lagrangian. Moreover we are in the non-perturbative regime 
of the strong interaction {E < Mj- ~ 1.7 GeV), so any perturbative treatment within QCD will 
not be a good approach to the problem. 

At very low energies {E < Mp, being Mp the mass of the p{770) resonance), the chiral symmetry 
of massless QCD rules the construction of an effective field theory that allows a perturbative 
expansion in momenta (p) and light quark masses (m), as (p^, m'^)/ A^, being A-^ ~ AttF ~ Mp 
the chiral symmetry breaking scale. rriT^ is the pion mass and F is its decay constant, with the 
normalization of F = 92.4 MeV. This theory is xPT, that drives the hadronization of QCD 
currents into the lightest multiplet of pseudoscalar mesons, including vr, K and rj particles. This 
framework was applied to the two- and three-pion tau decays in Ref. [2G] and it was checked that 
it provides a right description of the low-energy data though it was incapable of providing a good 
parametrization for the rest of the spectrum. This study put forward whatever structure due to 
resonance exchange that the form factors may have should match the chiral behaviour in the limit 
where xPT applies. As we brought forward before, the Kiihn-Santamaria models fail to fulfill that 
condition already at O(p^) in the chiral expansion [8, 21]. 

We shall extend to energies E > Mp, where the expansion parameter of is no longer 

valid. In fact, it is not known how to develop a dual effective theory of QCD in this region. 
However, there is a construct that has proven to be useful in this regime shedding light on the 
appropriate structure of the Lagrangian theory that we could use. This is yielded by the large- A'^c 
limit of SU{Nc) QCD [27] which introduces the inverse of the number of colours (three in the real 
world) as the parameter to build the expansion upon. In our context, it amounts to consider an 
spectrum of an infinite number of zero-width resonances that interact at tree level through local 
effective vertices. This frame, as we will see, can be used to establish a starting point in the study 
of the resonance region, and consequently in the hadronic decays of the tau lepton we are dealing 
with. The setting recalls the role of the Resonance Chiral Theory (RxT) [28, 29] that can be better 
understood in the light of the large- Ac limit of QCD as a theory of hadrons[30, 31]. 

The theory, although built upon symmetries guided by the large- Ac expansion and reproducing 
the chiral behaviour in the low-energy region, is still missing an ingredient of QCD. At higher 
energies {E 3> Mp), where the light-flavoured continuum is reached, perturbative QCD is the 
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appropriate tool to deal with the description of interactions, which is given in terms of partons. 
A well-known feature of form factors of QCD currents is their smooth behaviour at high energies 
[17, 18]. Then, it is plausible that the form factors match this behaviour above the energy region 
where the resonances lie. A complementary approach is given by the study of the operator product 
expansion {OPE) of Green functions of QCD currents that are order parameters of the chiral 
symmetry breakdown. Refs. [32-37] have evaluated these Green Functions within a resonance 
theory and proceeded to match the outcome to the leading term of the OPE at high transfers of 
momenta. As we commented it is crucial to take into account this high-energy information to settle 
a resonance Lagrangian that implements as many QCD features as possible, whence the described 
procedure will help to establish relations among some of the Lagrangian couplings, and eventually 
fix others [•■)7] . 

We will consider the SD description of the processes of t~ — )■ P~^Vr (P = vr, K) within the 
framework of RxT as introduced in this section and discussed in detail in the following one. These 
channels have not been observed yet, which is strange according to the most naive expectations of 
their decay rates. Clarifying this question is one of the main motivations of our study. 

The relative sign between the IB and SD dependent part motivated an addendum to Ref. [4]. 
This confusion was caused by the fact that they did not use a Lagrangian approach. In any such 
approach this should not be an issue. In order to facilitate any independent check, we define the 
convention we follow as the one used by PDG [38] in order to relate the external fields r^,^^ with 
the physical photon field 

r^, = i^, = -eQA^, (1) 

where e is the electric charge of the positron and the charge matrix of the three light flavour quarks 
\s Q = diag(|, — |, — |). Determining the relative sign between the structure independent {SI) and 
dependent contributions is an added interest of our computation. 

The SI part of the process has been discussed in [1]. We will compute the SD part using the 
Lagrangians in Sect. II. The kinematics and differential decay rates are discussed in Sect. Ill, as 
well as the general form of the matrix elements and the spectra. The structure dependent form 
factors for each mode are computed in Sects. IV and V. The QCD short-distance constraints are 
presented in Sect. VI and the phenomenological implications are detailed in Sect. VII. Finally, we 
summarize and discuss our results in Sect. VIII. 
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II. THEORETICAL FRAMEWORK 

The hadronization of the currents that rule semileptonic tau decays is driven by non-perturbative 
QCD. As mentioned in the Introduction, our methodology stands on the construction of an 
action, with the relevant degrees of freedom, led by the chiral symmetry and the known asymptotic 
behaviour of the form factors and Green functions driven by large Nq QCD. We will present here 
those pieces of the action that are relevant for the study of one meson radiative decays of the tau 
lepton. Hence we will need to include both even- and odd-intrinsic parity sectors. 

The large Nq expansion of SU{Nc) QCD implies that, in the Nq — )• oo limit, the study of 
Green functions of QCD currents and associated form factors can be carried out through the tree 
level diagrams of a Lagrangian theory that includes an infinite spectrum of strictly stable states 
[27] ^. Therefore the study of the resonance energy region can be performed by constructing such 
a Lagrangian theory. However, it is not known how to implement an infinite spectrum in a model- 
independent way. Moreover, it is well known from the phenomenology that the main role is always 
played by the lightest resonances. Accordingly it was suggested in Refs. [28, 29] that one can 
construct a suitable effective Lagrangian involving the lightest multiplets of resonances and the 
pseudo-Goldstone bosons (vr, K and ?/). This is indeed an appropriate tool to handle the hadronic 
decays of the tau lepton and the pion form factors [8-12, 40-42]. The guiding principle in the 
construction of such a Lagrangian is chiral symmetry. When resonances are integrated out from 
the theory, i.e. one tries to describe the energy region below such states {E <C Mp), the remaining 
setting is xPT, reviewed in Refs. [43, 44]. 

The very low-energy strong interaction in the light quark sector is known to be ruled by the 
SU{3)l (8> SU{3)r chiral symmetry of massless QCD implemented in xPT. The leading even- 
intrinsic-parity OijP') Lagrangian, which carries the information of the spontaneous symmetry 
breaking of the theory, is : 

C-% = ^{u,u^ + X+) , (2) 

where 

= i[u\df, - irp)u - u{dp - U^j,)u^] , 
X± = u^XU^ ± ux^u , X = '^Bq{s + ip) , (3) 

^ Since light resonances reach their on-shell peaks in the energy region spanned by the considered decays, the 
corresponding ofT-shell widths, which are energy-dependent, need to be implemented as we do following Refs. 
[12, 39]. 
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and (. . .) is short for the trace in the flavour space. The Goldstone octet of pseudoscalar fields 



is realized non-linearly into the unitary matrix in the flavour space 

= exp|-^$(x)| , (5) 

which under chiral rotations transforms as 

u{^) ^ 9Ru{^)h{9,^)^ = h{9,'P)u{ip) g[ , (6) 

with g = {gL,gR) G 5'[/(3)l ® 5J7(3)r and h{g,ip) G SU{3)v- External hermitian matrix fields 
''/i) ^/i) s and p promote the global 5'C/(3)l C5 S'f/(3)R symmetry to a local one. Thus, interactions 
with electroweak bosons can be accommodated through the vector v^^ = (r^ + i^)/2 and axial- 
vector = (r^ — l^)/2 fields. The scalar field s incorporates the explicit chiral symmetry breaking 
through the quark masses taking s = M + . . ., with M = diag(m„, nid, mg), where we will always 
work in the isospin limit in the present discussion, i.e. m„ = ma- Finally, at lowest order in the 
chiral expansion F is the pion decay constant and BqF^ = — (-i/i^), with the quark condensate of 
{tpip) = {uu) = (dd) = (ss). 

The leading action in the odd-intrinsic-parity sector arises at 0{p'^). This is given by the chiral 
anomaly [45] and explicitly stated by the Wess-Zumino-Witten (WZW) functional Zy[rzw[v,(i] 
that can be read in Ref. [43]. This contains all anomalous contributions to electromagnetic and 
semileptonic meson decays. For completeness, we give the relevant terms to r — ?• Pji'r below 

Cwzw = -^e^^afsi^^ W d-'r^ Uf + T.^^ ^ d^l^ + d-'l^ f - {L ^ R)) , (7) 

with U = u^, = U^d^U and = U^^'U^. 

It is well known [28, 46] that higher orders in the chiral expansion, i.e. even- intrinsic-parity 

(n) 

xPT '^ith n > 2, bring in the information of heavier degrees of freedom that have been integrated 
out, for instance resonance states. As our next step intends to include the latter explicitly, we 
will not consider higher orders in xPT in order to avoid double counting issues. In order to fulfill 
this procedure — at least, up to 0{p'^) in the even-intrinsic-parity sector — it is convenient to use 
the antisymmetric tensor representation for the J = 1 fields, as we comment below. Analogously, 
additional odd-intrinsic-parity amplitudes arise at 0{p^) in xPT, either from one-loop diagrams 



7 



using one vertex from the WZW action or from tree- level operators [47]. However we will assume 
that the latter are fully generated by resonance contributions [■'> I] and, therefore, will not be 
included in the following. 

The formulation of a Lagrangian theory that includes both the octet of Goldstone mesons and 
U (3) nonets of resonances is carried out through the construction of a phenomenological Lagrangian 
[48] where chiral and discrete symmetries determine the structure of the operators. Given the vector 
character of the Standard Model (SM) couplings of the hadron matrix elements in r decays, form 
factors for these processes are ruled by vector and axial-vector resonances. Notwithstanding those 
form factors are given, in the r — t- P^Vt decays, by a three-point Green function where other 
quantum numbers might play a role, namely scalar and pseudoscalar resonances [49]. However 
their contribution should be very small for r — t- P^fr- This statement is based on the following 
observations: the scalar resonances will be irrelevant at tree level to the considered process due to 
the discrete symmetry; about the pseudoscalar resonances, their contributions are suppressed due 
to first their heavy masses and also the fact that the relevant couplings involving the pseudoscalar 
resonances should be very tiny, since the decay of these states to P7 has not been reported yet. 
Thus in our description we include J = 1 resonances only, and this is done by considering a nonet 
of fields [28] : 

1 * 

R = ^Y.^i't>R,^ , (8) 

where R = V,A, stand for the vector and axial- vector resonance states. Under the SU{2>) l®SU{^)r 
chiral group, R transforms as : 

R ^ h{g,^)Rh{g,^)'^ . (9) 

The fiavour structure of the resonances is analogous to that of the Goldstone bosons in Eq. (4). 
One can also introduce the covariant derivative 

V^X = d^X + [T^,X] , (10) 
= ]^[v){d^-ir^)u + u{d^-U^)u^] , 

acting on any object X that transforms as R in Eq. (9), like and x±- The kinetic terms for the 
spin-one resonances in the Lagrangian read [28] : 

1 A/f^ 

jOL = -^{'^'Rx^VuR''n + ^{R,.uR^n , R = V,A, (11) 



8 



My, Ma being the masses of the nonets of vector and axial-vector resonances in the chiral and 
laige-Nc hmits, respectively. Notice that we describe the resonance fields through the antisym- 
metric tensor representation. With this description one is able to collect, upon integration of the 
resonances, the bulk of the low-energy couplings at 0{p^) in xPT without the inclusion of addi- 
tional local terms [oT] . So it is necessary to use this representation if one does not include ^^^pt 
in the Lagrangian. Though analogous studies at higher chiral orders have not been carried out, 

(n) 

we will assume that no ^^p-p with n = 4,6,... in the even-intrinsic-parity and n = 6,8,... in the 
odd-intrinsic-parity sectors need to be included in the theory. 

The construction of the interaction terms involving resonance and Goldstone fields is driven by 
chiral and discrete symmetries with a generic structure given by : 

Oi ^ {RiR2...RjX^''H^)), (12) 

where X^^\v^) is a chiral tensor that includes only Goldstone and auxiliary fields. It transforms 
like R in Eq. (9) and has chiral counting n in the frame of xP'^- This counting is relevant in the 
setting of the theory because, though the resonance theory itself has no perturbative expansion, 
higher values of n may originate violations of the proper asymptotic behaviour of form factors 
or Green functions. As a guide we will include at least those operators that, contributing to our 
processes, are leading when integrating out the resonances. In addition we do not include operators 
with higher-order chiral tensors, that would violate the QCD asymptotic behaviour unless 

their couplings are severely fine tuned to ensure the needed cancellations of large momenta. In the 
odd-intrinsic-parity sector, which contributes to the vector form factor, this amounts to include 
all (-Rx^^^> and (RRx^^^) terms. In the even- intrinsic-parity sector, contributing to the axial- 
vector form factors, these are the terms {Rx^'^'^)- However previous analyses of the axial- vector 
contributions [8, 11, 12, 35] show the relevant role of the {RRx^'^^) terms that, accordingly, are also 
considered here. 

We also assume exact SU (3) symmetry in the construction of the interacting terms, i.e. at the 
level of couplings. Deviations from exact symmetry in hadronic tau decays have been considered in 
Ref. [50] . However we do not include SU{3) breaking couplings because we are neither considering 
next-to- leading order corrections in the 1/Nc expansion. 

The lowest order interaction operators, linear in the resonance fields, have the structure 
{Rx'-'^H^)). There are no odd- intrinsic-parity terms of this form. The even-intrinsic-parity La- 
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grangian includes three coupling constants [28] : 



= {^{^uD , (13) 

where f!^ = uF'^^u^ -iiu^ F^^u and F'^j^ are the field strength tensors associated with the external 
right- and left-handed auxiliary fields. All couplings Fy, Gy and Fa are real. 

The leading odd-intrinsic-parity operators, linear in the resonance fields, have the form 
(Rx^^^if))- We will need those pieces that generate the vertex with one vector resonance, a 
vector current and one pseudoscalar. The minimal Lagrangian with these features is : 

where Cj are real dimensionless couplings, and the VJP operators read [34] 

O^jP = ie;..p.(^^n£",X+]> , 

Ol, = e,..p.({V-y^^/f}n,) . (15) 

Notice that we do not include analogous pieces with an axial-vector resonance, that would con- 
tribute to the hadronization of the axial- vector current. This has been thoroughly studied in 
Ref. [n] (see also Ref. [12]) in the description of the r — t- vrvrvrz/T- process and it is shown that no 
(Ax^'^^f)) operators are needed to describe its hadronization. Therefore those operators are not 
included in our minimal description of the relevant form factors. 

In order to study tau decay processes with a pseudoscalar meson and a photon in the final state 
one also has to consider non-linear terms in the resonance fields. Indeed the hadron final state in 
T — )• Pji'r decays can be driven by vertices involving two resonances and a pseudoscalar meson. 
The structure of the operators that give those vertices is {RiR2X^'^\^))i and has been worked out 
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before [8, 34]. They include both even- and odd-intrinsic-parity terms : 

5 4 

cr = E^^^vAP + E'^^^vvp, (16) 

i=l i=l 

where Aj, and di are unknown real dimensionless couplings. The operators O^j^p are given by : 
1/ VAP terms 



(17) 



'-^VAP 


= {[V^\A^,\x-) , 






03 

*-^VAP 
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with h^j_y = V^Uj, + VuUfj., and 
2/ VVP terms 

We emphasize that is a complete basis for constructing vertices with only one pseudoscalar 
meson; for a larger number of pseudoscalars additional operators might be added. As we are only 
interested in tree-level diagrams, the equation of motion arising from ^^px ™ (^) been 
used in CJ and to eliminate superfluous operators. 
Hence our theory is given by the Lagrangian : 

^RxT = ^xPT + ^WZW + -^kin + C2 + C2 + Cj + . (19) 

It is important to point out that the resonance theory constructed above is not a theory of QCD for 
arbitrary values of the couplings in the interaction terms. As we will see later on, these constants 
can be constrained by imposing well accepted dynamical properties of the underlying theory. 
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III. MATRIX ELEMENT DECOMPOSITION, KINEMATICS AND DECAY RATE 

The process we are going to compute is T~{pr) — )• Vt{<1)P~ {p)l{k) ■ The kinematics of this 
decay is equivalent to that of the radiative pion decay [51]. We will use t := {p^ — q)^ = {k + 
In complete analogy to the case of the radiative pion decay [52], the matrix element for the decay 
of T~ — )• P~^Vt- can be written as the sum of four contributions: 

M[T~{pr)^yr{q)P~{p)l{k)\ =Mib^+Mibp+Mv+Ma, (20) 

with ^ 

iMiBr = GFV^'j,^jeFpp^e,{k)L^'\ 
iMv = iGFV'c'KM^Pv{t)e^,upae''{k)kPp'' L^' , 

iMA = Gf V'Skm e FA{t) e'ik) [{t - m|) g^, - 2p, k^] L^' , (21) 

where is the polarization vector of the photon. Fv{t) and FA{t) are the so-called structure 
dependent form factors. Finally and L'^'^ are lepton currents defined by 

= U,^{q)Yil-75)Ur{Pr): 

L^" = u.r{<l)l'{l - 75) (/_"^^;^,"_^^, 7-n.(p.) . (22) 

The notation introduced for the amplitudes describes the four kinds of contributions: Mis-r is 
the bremsstrahlung off the tau, (Figure 1(a)); MiSp is the sum of the bremsstrahlung off the P- 
meson (Figure 1(b)), and the seagull diagram (Figure 1(c)); My is the structure dependent vector 
contribution (Figure 1(d)) and Ma the structure dependent axial- vector contribution (Figure 1(e)). 
Our ignorance of the exact mechanism of hadronization is parametrized in terms of the two form 
factors FA{t) and -Fy(t). In fact, these form factors are the same functions of the momentum 
transfer t as those in the radiative pion decay, the only difference being that t now varies from 
up to rather than just up to m^. 

The two matrix elements Mi Br ^^^^ Mjbp are not separately gauge invariant, but their sum, 
ie. the (total) matrix element for internal bremsstrahlung IB 

MiB = MiBr+MiBp, (23) 



^ Notice that i and minus factors differ with respect to Ref. [1] (DF). Moreover, our form factors have dimension of 
inverse mass while theirs are dimensionless. In their work, the factor of {y/2m^)~^ in the form factors is compen- 
sated by defining tlie sum over polarizations of the matrix element squared with an extra 2m^ factor. This should 
be taken into account to compare formulae in both works using that = ^/2m^Fv(^)°"^ Fa^'''' = 
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(a) -y Vr (b) i/r -y (c) Vr . 7 

■ /a"a . £: ■ /a"a . ^ P- C_ 




(d) 1/ 7 _ (e) fr, 7 

jv: 





FIG. 1: Fcynman diagrams for the different kinds of contributions to the radiative decays of the tau including 
one meson, as explained in the main text. The dot indicates the hadronization of the QCD currents. The 
solid square represents the SD contribution mediated by the vector current and the solid triangle the SD 
contribution via the axial- vector current. 

is indeed gauge invariant, as Aiy and Ma are. We also define the (total) structure dependent 
radiation SD hy 

Msd = Mv+Ma- (24) 
The spinor structure can be rearranged to give 



, , n , '^APt), (25) 

Pr • k p ■ k Zpr • k J 

iMsD = Gf F.^"^*' e {ie^,p„L^'e''kPp''Fv{t) + u,M){^ + 75) [{t - ml)^ - 2(e • p)^] u{pr)FA{t)] . 
The square of the matrix element is then given by 



lA^P = \MiB? + l^eiMiBM^so) + \Msd? , (26) 

where the bar denotes summing over the photon polarization, the neutrino spin and averaging over 
the tau spin. 

We follow Ref. [4] and divide the decay rate as follows: the internal bremsstrahlung part Tjb 



arising from jA^/^p, the structure dependent part Tsd coming from I^WsdI ) and the interference 



part Tint stemming from 2^e{A4 jb-Msd)- Furthermore Tsd is subdivided into the vector-vector 
(Fyy), the axial-vector-axial-vector (Taa) and the vector-axial-vector interference term (Tva)- 
Similarly F/atj' gets split into the internal bremsstrahlung- vector interference (Tjb-v) ^-nd the 
internal bremsstrahlung-axial-vector interference {Tjb~a) parts. Thus, one has 

Tall = Tjb + Tsd + Tint j 
Tsd = Tvv + Tva + Taa , 

Tint = Tib-v + Tib-a ■ (27) 
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It is convenient to use the dimensionless variables x and y to proceed, as used in Ref.[l] and 
references therein: 

In the tau rest frame x and y are the energies and of the photon and the pion, respectively, 
expressed in units of Mt-/2: 

E^ = ^x, E^ = ^y. (29) 

Eq. (29) sets the scale for the photons to be considered as "hard" or "soft". This means that the 
formulae for internal bremsstrahlung should be similar for radiative tau and pion decay, once they 
are expressed in terms of x and y, as it is the case, albeit photons of comparable softness will have 
very different energies in both cases. 

The kinematical boundaries for x and y are given by 

0<x<l-r|,, l-x + ^<y<l + rl, (30) 

where 

2 _ / 'rnp\ 0.006^ 1 fo-^s 

\ Tfj ~o.077< 1 > (31) 

where the upper figure corresponds to P = tt and the lower one to P = i^. It is also useful to note 
that 

p.k=^{x + y-l-rl), t:={pr-q)^ = ik+pf = M^{x + y-l). (32) 
The differential decay rate is given by [53] 

^<-' - = '''' 

where the bar over the matrix element denotes summing over the photon polarization, the neutrino 
spin and averaging over the tau spin. Choice of the tau rest frame, integration over the neutrino 
momentum, p, and the remaining angles and introduction of x and y yield 

(Pr rrir 



dx dy 2567r3 
The integration over y yields the photon spectrum 



^l-MP- (34) 



dr r^^^'p d^r 

^ = / ^ ■ (35) 

dx Ji_^+lp_ dxdy 

1 — x 
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A low-energy cut must be introduced for the photon energy because of the infrared divergence of 
the internal bremsstrahlung. By requiring x > xq one obtains the integrated decay rate 

/■I -4 dr 

r(xo) = T{Eo) = J dx— , (36) 

that does depend on the photon energy cut-off {Eq = ^xq). Instead of x and y one can also use 
X and z, where z is the scaled momentum transfer squared: 

^ = ]j^ = ^ + y- l' (37) 

whose kinematical boundaries are 

z-rl<x<l-^, rl< z<l. (38) 

Integration of over x yields the spectrum in z, i. e. the spectrum in the invariant mass of the 
meson-photon system: 

dr, , dr / ^\ , d^r , 

The integrated rate for events with t > to is then given by 



r(zo) = r (VTo) = f 

J Zn 



dz'^iz). (40) 



We note that zq is both an infrared and a collinear cut-off. 

In terms of the quantities defined in Eq. (27) the differential decay rate is 



d r/B a , 



fjB {x,y,r%) \^~^''f- , (41) 



dxdy 2vr^^-^"^'"'--^ (i_^2)^ 

[\Fv{t)\^fvv {x,y,rl)+meiFv{t)FXit))fvA {x,y,rl) + 



d^TsD _ a 



dxdy Svr Fp 

4\FA{t)\'fAAix,y,rl)]^-^_^^^^^ , 

= £^ [fiB^v {x,y,r%) ^eiFvit)) + 2fjB-A {x,y,rl) ^FAit))] ^^^'^'"/^ ' 

with a = e^/(47r) standing for the fine structure constant and 

[rj,{x + 2) - 2rj,{x + y) + {x + y - 1) {2 - 3x + x'^ + xy)] {rj, - y + l) 



fiB {x, y, rp) 



rp — X — y + I] X 



fvv{x,y,rp) = -[rj,{x + y) + 2rp{l-y){x + y) + {x + y-l){-x + x'^ -y + y^)], 

fAA {x, y, rp) = fvv {x, y, rp) , 

fvA {x, y, rp) = - [rp{x + y) + {I - x - y){y - x)] {rj, - x - y + l) , 

/ 2^ {rl-x-y + 1) {rl-y + l) 

fiB-v [x,y,rp) = ^ -, 

2^ [rj,-2rl{x + y) + {l-x + y){x + y-l)]{rl-y+l) 

fiB-A{x,y,rp) = ^ —r- . (42) 

[Tp — X — y + 1) X 
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The radiative decay rate has been expressed in term of the width of the non-radiative decay 

r^2 12 772 



We finish this section by presenting the analytical expressions for the invariant mass spectrum: 
^^^^ - ^ [,.4(1 + 2rl [z - z') - 4z + 5z^ - z' 

1 r 



dz 2tt 



+ [r%z + 2rj,z - 2z - 2z^ + z^) Inz] 



^2_^2w. __5^2 ' 



"P^ (1-4) 

dz 247rF2 ^2 l^l^WI (i_^2^)2 

dz ~ evrFl z2 I^^WI (i_^2^)2 ' 

dr/p„y _ a M'^{z-rlf{l-z + z\nz) T,-^^^p- 



^e{Fv{t))- 



dz 2-K Fp z (^i _ 



dT,s-A _ -^ir^rl{l-z)-l-z + 2z- 



dz IT Fp 



+ ^rlz -2z- z') Inz] ^^jRe(FA(t)) 7'^-;^, . (44) 

z (1 - r%) 

The interference terms IB — V and IB — A are now finite in the limit z — )• rp, which proves that 
their infrared divergences are integrable. 

Although the above formulae have been noted in Ref.[4], we independently calculate them and 
explicitly give them here for completeness. Moreover we would like to point out that due to the 
fact that our definitions of the form- factors Fv{t) and F^^t) differ from the ones given in Ref.[4], 
as we have mentioned before, there are some subtle differences in the above formulae between ours 
and theirs. 

IV. STRUCTURE DEPENDENT FORM FACTORS IN r ^ tt 71/^ 

The Feynman diagrams, which are relevant to the vector current contributions to the SD part 
of the T~ — )■ ■K~^Vr processes are given in Figure 2. The analytical result is found to be 

iMsT,y=iGFVudeuuM)r{^-l5)ur{s)e^,o.pe^{k)k''p^F^{t), (45) 



We disagree with Ref. [i] on the signs of fvA and fiB-v in Eq.(42) and — ^ in Eq.(44), even after taking 
into account the minus sign difference in the definition of the IB part. 
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where the vector form- factor Fy[t) is 

Nr. 

mt) -- 



Nc ^ 2V2Fv 



(C2 - Ci - C^)t + (C5 - Ci - C2 - 8C3)m^ 



cos 



+ 



1 - V2tge] + 



sm 



M2 



1 + \/2cotg(9 



(Ci - C2 - C5 + 2c6)t + (C5 - Ci - C2 - 8c3)m^ 



4F, 



cos^^ 



Mi 



1 - V2t^e] + 



sm 



M.2 



1 + V2cotg(9 



Here we have defined Dnit) as 



(46) 



(47) 



Ml - t - iMRTR{t) ■ 

Vii[t) stands for the decay width of the resonance R, which will be included following Refs. [12, 39]. 
For completeness, we write the explicit expressions of the off-shell widths in Appendix A. 

We will assume the ideal mixing case for the vector resonances a; and in any numerical 
application: 



oji = cosO Lo — sm.9 



UJ 



2 1 
V 3 V 3 



(48) 



WW ®VA/W\A/ (^VA/WW (^VA/WVWW 

FIG. 2: Vector current contributions to r~ — TT~^i'r- 

The Feynman diagrams related to the axial- vector current contribution to the SD part are given 
in Figure 3. The corresponding result is 

iMsD^=GFVudeUuAQ)l''i^-75)uris)e''ik) [{t - ml)g^, - 2k^p,] F^it) , (49) 
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where the axial- vector form- factor F'^{t) is 



2Gv\ Fl 



2F^Mj \} Fv ) 2F^ 



V2FaFv 
F^Mj 



D,,{t) -X"t + Xoml], 



(50) 



where we have used the notation 



V2X0 
V2X" 



-4Ai - A2 - y - A5 , 
A2-y-A5, 



for the relevant combinations of the couplings in £.2"^^ ^ Eq. (17). 



(51) 



p" 7 



p 



(gVAA#WV ^VA/WVWW 



FIG. 3: Axial-vector current contributions to r — > tt 71^1- 



V. STRUCTURE DEPENDENT FORM FACTORS IN r 



K JVr 



Although one can read the following observation from Eq.(21), let us emphasize that the model 
independent part M.ib^+k same as in the pion case by replacing the pion decay constant 

with the kaon decay constant Fx- A brief explanation about this replacement is in order. The 
difference of F^^ and Fx is generated by the low-energy constants and the chiral loops in xPF [15], 
while in the large Nc limit of RxT this difference is due to the scalar resonances in an implicit way. 
Due to the scalar tadpole, one can always attach a scalar resonance to any of the pseudo-Goldstone 
boson fields, which will cause its wave function renormalization. A convenient way to count this 
effect is to make the scalar field redefinition before the explicit computation to eliminate the scalar 
tadpole effects. In the latter method, one can easily get the difference of and Fx- For details, 
see Ref. [54] and references therein. For the model dependent parts, the simple replacements are 
not applicable and one needs to work out the corresponding form factors explicitly. 

The vector current contributions to the SD part of the t~ — >• K~^Ur process are given in Figure 
4. The analytical result is found to be 



iMsjiy = iGf Vus eUuM) T^'Cl - 75) Ur{s)e^^o.p e'ik) k'^p^ F^ {t) 



(52) 
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where the vector form- factor Fy(t) is 



Nc V2F^ 



+ 



24tt^Fk FkMv 



1 - 2V2cotg(9 



M2 3M2 



3FkMv ^ ' 



(C2 - Ci - C5)t + (C5 - Ci - C2 - 8C3)m|^ 

1 + 2v^tgi 



3Mi 



+ 



1 sin^ 



M2 3M2 



Ci - C2 - C5 + 2c6)t + (C5 - Ci - C2 - 8c3)m|- 
dsi + (di + 8d2 - d3)mj^ 



2^2 



1 - 2V2cotg6l 



3Mi 



1 + 2^/2tg( 



(53) 



K- 



®VAAy 



A'- 



K*- 



K- 



K- 



K*- 0) 



FIG. 4: Vector current contributions to r K ji/r- 



The axial-vector current contributions to the SD part are given in Figure 5. The corresponding 
analytical result is 

iMsD^ = Gf Vus eu,M) 7^(1 - Ts) Ur{s) e'ik) [{t - m\)g^, - 2k^p,] F^{t) , (54) 
where the axial- vector form-factor F^{t) is 



sin2( 



+ 



M2, 



2Fk 



+ 



FaFv 
V2Fk 



1 



+ 



cos20 sin2^ 
+ 



M2 Ml M2 



A"t -|- \Qm\ 



(55) 



We have used the notations of Kih and Kn for the physical states i('i(1400) and i('i(1270), 
respectively, and the mixing angle 6a is defined in Eq.(56) as we explain in the following. 
The KiA state is related to the physical states iCi(1270), -ftri(1400) through: 



KiA = cosOa i^i(1400) -Fsin6lA i^i(1270) . 



(56) 



The nature of Ki[12'l{)) and /^i(1400) is not clear yet. It has been proposed in Ref. [55] that they 
result from the mixing of the states Kia and Kib , where Kia denotes the strange partner of the 
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axial vector resonance ai with J^^ = and Kib is the corresponding strange partner of the 
axial vector resonance hi with J^^ = l'^ . However in this work, we will not include the nonet of 
axial vector resonances with J^*^ = 1"^ V'^']- As argued in Ref. [55], the contributions from these 
resonances to tau decays are proportional to SU (3) symmetry breaking effects. Moreover, as one 
can see later, we will always assume SU (3) symmetry for both vector and axial-vector resonances 
in deriving the T-matrix. For the pseudo-Goldstone bosons, physical masses will arise through 
the chiral symmetry breaking mechanism in the same way as it happens in QCD. For the vector 
resonances, the experimental values will be taken into account in the kinematics, while in the case 
of the axial-vector mesons, we will take the determination of the ai mass from Ref. [12] and the 
masses of the Ki resonances from PDG [38]. 



K- 

/ 7 



K- 



K- 



K- 



FIG. 5: Axial- vector current contributions to t — > X ji^r- 



VI. CONSTRAINTS FROM QCD ASYMPTOTIC BEHAVIOUR 

In this part, we will exploit the asymptotic results of the form factors from perturbative QCD 
to constrain the resonance couplings. When discussing the high-energy constraints, we will work 
both in chiral and SU (3) limits, which indicates we will not distinguish the form factors with pion 
and kaon, that are identical in this case ^. 

For the vector form factor, the asymptotic result of perturbative QCD has been derived in 
Ref. [18, 57] 

FP{t^-oo) = j, (57) 

where F is the pion decay constant in the chiral limit. From the above asymptotic behaviour, we 
find three constraints on the resonance couplings 

Cl - C2 + C5 = , (58) 



* The results of this procedure are independent of taking the chiral limit [34]. 
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C-c+c- 2C6 = 3,^,^^ + , (59) 

C2-C1 + C5-2C6 = ^^r;^ + ^^3, (60) 

where the constraints in Eqs.(58), (59) and (60) are derived from 0{t^), 0{t^) and 0{t~^), respec- 
tively. Combining the above three constraints, we have 

NcMy , Fy ^ 

where the constraint of Eq.(62) has already been noticed in [4, 18, 57]. 

It is worthy to point out different results for the asymptotic behavior of the vector form factor 
F-,^^[t) have also been noted in different frameworks, such as the ones given in Refs. [58, 59]. In 
Refs. [18, 57], the result was obtained in the parton picture and unavoidably the parton distribution 
function for the pion has to be imposed to give the final predictions. OPE technique was exploited 
to obtain its prediction in [5S], which led to the conclusion that potentially large QCD corrections 
could exist. In [59], the form factor was discussed by using Bjorken- Johnson-Low theorem [60]. 
Variant methods have also been used to analyze this form factor: the corrections from the transverse 
momentum of the parton were addressed in Ref. [61]; a QCD sum rule method was applied to 
derive the asymptotic behavior in Ref. [62]. In the present discussion, we focus our attention on 
Refs. [58, 59], while the study for the other results can be done analogously. Although different 
results agree with the same leading power of the square momentum for large t, behaving as 1/t, 
they predict different coefficients, such as 

F^(t^-oo) = — , (63) 

from Ref. [58] and 



FP{t ^ -oo) = I . (64) 

from Ref. [59]. 

By doing the same analyses as we have done by using the Lepage-Brodsky result in Eq.(57) 
to constrain the resonance couplings, we can straightforwardly get the constraints from the short- 
distance behaviors given in Eqs.(63)-(64). Apparently the matching results given in Eqs.(58)-(59) 
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will stay the same, since they are derived from 0{t^) and 0{t^). Comparing with the result of 
Eq.(60) from the matching of 0{t~^) by using the coefficient in Eq.(57), the corresponding results 
by using Eqs.(63) and (64) are respectively 

c,-c, + C5-2c6 = ^7|^ + ig^rf3, (65) 

c.-c, + C5-2ce = ^7|^ + 1^^^3, (66) 

which lead to the following results, in order, by combining Eqs.(58)-(59) 

,67, 

The formulae displayed in Eq.(62), Eq.(67) and Eq.(68) provide a simple way to discriminate 
between different asymptotic behaviors. The chiral limit values for the pion decay constant F 
and the mass of the lowest vector multiplet My have been thoughtfully studied at the leading 
order of l/Nc in Ref. [03], which predicts F = 90.8 MeV and My = 764.3 MeV. The different 
results shown in Eq.(62), Eq.(67) and Eq.(68) from matching different short-distance behaviors of 
Refs. [18, 58, 59] deviate from the phenomenology study at the level of 5%, 16%, 64% respectively 
^, which implies that the short-distance behavior in Eq.(57) is more reasonable than the ones in 
Eqs.(63)-(64). Hence for the matching result at the order of 1/t, we will use the one in Eq.(62) 
throughout the following discussion. However we stress the inclusion of extra multiplets of vector 
resonances or the sub-leading corrections in 1/Nc may alter the current conclusion. 

The high-energy constraints on the resonance couplings Cj and di have been studied in different 
processes. The OPE analysis of the VVP Green Function gives [34] 



NcM^ , F' 



The constraint from r — ?• (VP) study leads to 



\/2My 

if one neglects the heavier vector resonance multiplet [10]. 



C5 - C6 = ds , (71) 



^ The results are mildly changed when estimating the values of F and My by and Mp respectively, as expected. 
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The results from the analysis of r — )■ (KKtt) Uj- are [64] 

NcMvFv 



C5 - C6 



192^/27r2F2 ' 
1927r2F2 ■ 

It is easy to check that the results of Eqs.(71) and (72) are consistent. Combining Eqs.(61) and 



_ NcMl 



(71) leads to 



C5 - C6 



64V27r2Fv/ ' 



where the constraint of C5 — cg is consistent with the result from the OPE analysis of the VVP 
Green Function [ : ], while the result of ^3 is not. 

By demanding the consistency of the constraints derived from the processes of — )• P^jv^- 
and T~ —7- {VP)^i'r given in Eq.(73) and the results from t~ — > {KKiT)~Vr given in Eq.(72), we 
get the following constraint 

Fv = VSF . (74) 

If one combines the high-energy constraint from the two pion vector form factor [37] 

FvGv = , (75) 

and the result of Eq.(74) we get here, the modified Kawarabayashi-Suzuki-Riazuddin-Fayyazuddin 
(KSRF) relation [65, 66] is derived 

F = VSGv , (76) 

which is also obtained in the partial wave dispersion relation analysis of vrvr scattering by properly 
including the contributions from the crossed channels [67] . 

Although the branching ratios for the modes of r — t- Pji'r which we are discussing could 
be higher than for some modes that have been already detected, they have not been observed 
yet. Lacking of experimental data, we will make some theoretically and phenomenologically based 
assumptions in order to present our predictions for the spectra and branching ratios. 

Taking into account the previous relations one would have Fy{t) in terms of ci + C2 + 8C3 — C5 
and di + 8^2 — d^- For the first combination, ci + C2 + 8C3 — C5 = Ci + 4c3 (ci — C2 + C5 = has 
been used), the prediction for c\ + 4c3 in [34] yields ci + C2 + 8C3 — C5 = 0. In Ref. [34] the other 
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relevant combination of couplings is also restricted: di + 8d2 — ds = -^bt- In Fy{t), C4 appears 
in addition. There is a phenomenological determination of this coupling in the study of the KKn 
decay modes of the r [11]: C4 = —0.07 it 0.01. 

Turning now to the axial-vector form factor, it still depends on four couplings in both channels: 
Fa, Ma, A" and Aq- If one invokes the once subtracted dispersion relation for the axial vector 
form factor, as done in Ref. [4], one can not get any constraints on the resonance couplings from 
the axial vector form factors given in Eqs.(50) and (55). In fact by demanding the form factor to 
satisfy the unsubtracted dispersion relation, which guarantees a better high-energy limit, we can 
get the following constraint 

which has been already noted in [-35]. 

In order to constrain the free parameters as much as possible, we decide to exploit the constraints 
from the Weinberg sum rules (WSR) [Gn]: F^ - F^ = F^ and M^F^ - M\F\ = 0, yielding 

Fa = 2F2 , Ma = -^. (78) 

For the axial vector resonance coupling Aq, we use the result from Ref. [12, 35] 

Ao = . (79) 

To conclude this section, we summarize the previous discussion on the high-energy constraints 

Fv = V3F, Gv = ^, Fa = V2F, My = V_ , Ma 



In the above results, we have discarded the constraint in Eq.(70), which is the only inconsistent 
result with the others. 



VII. PHENOMENOLOGICAL DISCUSSION 



Apart from the parameters we mentioned in the last section, there is still one free coupling 6a, 
which describes the mixing of the strange axial vector resonances in Eq.(56). The value of 9a has 
already been determined in the literature [10, 55, 69]. We recapitulate the main results in the 
following. 
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In Ref. [55], it has given Oa ~ 33°. In Ref. [10], \eA\ ~ 58.1° is determined through the 
considered decays r~ — )• {VP)~Vr. In Ref. [G9], the study of r — )• Kiv^- gives \6a\ =530 as the 
two possible solutions. The decay D — t- Kitt allows to conclude that 9a must be negative and 
it is pointed out that the observation of — )■ K^tt'^ with a branching ratio ~ 5 • 10~^ would 
imply 9a ~ —58°. However, a later analysis in Ref. [70] finds that the current measurement of 

—7- i^j^(1400)7r"'" [38] favors a mixing angle of —37° over —58°. In this respect, the relation 



r ( J/^' ^ K?(1400)K 



(81) 



would be very useful to get 9a, once these modes are detected. In the following discussion, we will 
show the results using both \9a\ = 37° and 58°. The other inputs are given in the Appendix.B. 



A. Results only with WZW contribution 

As it was stated before it is strange that the decay modes r — t- P^Vr have not been detected 
so far. The most naive and completely model independent estimate would be to just include the 
IB part and the WZW contribution to the VV part, as the latter is completely fixed by QCD. 
We know that doing this way we are losing the contribution of vector and axial- vector resonances, 
that should be important in the high-x region. However, even doing so one is able to find that the 
radiative decay t~ — t- tt~ji^t- has a decay probability larger than the mode t~ — t- K~ K~ Vt ^. 
For a reasonably low cut-off on the photon energy this conclusion holds for the t~ — )■ K"^Vt- as 
well. 

Before seeing this, we will discuss briefly the meaning of cutting on the photon energy As it is 
well known [71, 72] the IR divergences due to the vanishing photon mass cancel when considering 
at the same time the non-radiative and the radiative decays with one photon ^. In practice, this 
translates into mathematical language the physical notion that the detectors have a limited angular 
resolution that defines a threshold detection angle for photons. If one considers a photon emitted 
with a smaller angle it should be counted together with the non-radiative decay as it is effectively 
measured in this way. The sum is of course an IR safe observable. The splitting depends on the 
particular characteristics of the experimental setting. Obviously, the branching fraction for the 
radiative decay depends on this cut-off energy. We will consider here the case E^thr = 50 MeV, 

r (r" -> K+R-R-v-r) = 3.579(66) • 10~" GeV. 
^ A cut on the photon energy was introduced in Sect. III. 

* In general, the IR divergences of the n-photon decay are canceled by those in the n + 1-photon process. 
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that corresponds to x = 0.0565. In order to illustrate the dependence on this variable, we will also 
show the extremely conservative case of E^thr = 400 MeV (x = 0.45). In the first case we obtain 
r (r~ — )• 7r~7i/^) = 3.182 x 10~^^ GeV, and in the second one we are still above the 3K decay width, 
r (r~ — )• TT^^Ur) = 3.615 X 10^^^ GeV. In Figure 6 we plot the photon spectrum of r" — > ir^^Ur- 
Proceeding analogously for the decay with a K" , we find: T (r~ — t- K~'^Vr) = 6.002 x 10~^^ GeV 
for E^thr = 50 MeV, and T (r" ^ K'-fUr) = 4.589 x 10"^^ GeV for E^thv = 400 MeV. The 
photon spectrum of t~ — )• K~^Vr can be seen in Figure 7. For any reasonable cut on these 
modes should have already been detected by the i?-factories. 




FIG. 6: Differential decay width (in units of Tt^t^i,^ ) of the process — > ti^^Vt including only the model 
independent contributions as a function of x . For the form factors only the WZW term is considered for 
this estimate, where the axial vector contribution is absent. The right plot is the close-up of the left one in 
region of .t > 0.5. 



Already at this level of the phenomenological analysis, the question of the accuracy on the 
detection of soft photons at i?-factories [73] arises ^. An error larger than expected (here and 
in some undetected particle interpreted as missing energy, in addition to a gaussian treatment of 
systematic errors) could enlarge the uncertainty claimed on the measurement of B~ — t- t~z^t- [38] 
when combining the Belle [75] and BaBar measurements [76, 77] taking it closer to the Standard 
Model expectations. 



See however, Ref.[74] 
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FIG. 7: Differential decay width (in units of Fr-j-Kiy^ ) of tlie process K^ji^r including only the 

model independent contributions as a function of x . For the form factors only the WZW term is considered 
for this estimate, where the axial vector contribution is absent. The right plot is the close-up of the left one 
in region of a; > 0.5. 

B. Results with the resonance contributions for the pion channel 

Next we include also the model-dependent contributions. Since in the kaon channel there are 
uncertainties associated with the off-shell width of the strange axial-vector resonance and the 
mixing of the corresponding light and heavy states, we will present first the pion channel where 
there are not any uncertainties of these types and everything is fixed. 

In Fig. 8 the resulting photon spectrum of the process t~ — ?• tt^^Uj- is displayed. In order to 
display clearly how the different parts contribute to the spectrum, we have given the close-up of 
the spectrum for the high x region and also shown the separate contributions in Fig. 8. For "soft" 
photons {xq < 0.3) the internal bremsstrahlung dominates completely. One should note that for 
very soft photons the multi-photon production rate becomes important, thus making that our 
0{a) results are not reliable too close to the IR divergence x = 0. We qualitatively agree with 
the results in DF papers, for the same order of a to the three significant figures shown in Ref. [4]. 

The spectrum is significantly enhanced by SD contributions for hard photons (xq ^ 0.4), as we 
can see in Fig. 8. From Fig. 8 one can also see that the vector current contribution mediated by the 
vector resonances dominates the SD part. The interference terms between the bremsstrahlung and 
the SD parts are also shown in Fig. 8. If we compare the predicted curves with those in Ref. [4] we 
see that the qualitative behaviour is similar: the IB contribution dominates up to x ~ 0.75. For 
larger photon energies, the SD part is predominantly due to the VV contribution and overcomes 
the SI part. We confirm the peak and shoulder structure shown at x ~ 1 in the interference 
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contribution, which is essentiahy due to IB — V term, and also the V — A term, that is in any case 
tiny. 
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FIG. 8: Differential decay width ( in units oiTr^TTi'^ ) of the process r~ — > ■K~ji'r including all contributions 
as a function of x . The top right plot is the close-up of the top left one in region of a; > 0.5. The bottom left 
and right plots are to display the compositions of the interference and structure dependent parts, respectively. 



While the integration over the IB part needs an IR cut-off, the SD part does not. We have 
performed the integration over the complete phase space, yielding 

Tyy = 0.99 • 10-3 , TvA^O, T^a = 0.15 • IQ-^ ^ = 1.14 • 10"^ (82) 

Our number for Tsd lies between the results for the monopole and tripole parametrizations in Ref. 
[4]. However, they get a smaller(Iarger) yy(AA) contribution than we do by ~ 20%(~ 200%). 
This last discrepancy is due to the off-shell ai width they use. In fact, if we use the constant 
width approximation we get a number very close to theirs for the AA contribution. With the 
understanding of the ai width in the r — )• Siri'r observables [12], we can say that their (relatively) 



Here and in what follows, all contributions to the partial decay width are given in units of the non-radiative decay. 
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FIG. 9: Differential decay width ( in units of Fr-i-Tri/T ) of the process r n 71^1- including all contributions 
as a function of \/t. The top right plot is the close-up of the top left one in region 

of V< > 1-0. The middle 

left and right plots are to display the compositions of the interference and structure dependent parts, 
respectively. The bottom plot is the close up of the middle right one in the region of \/i > 0.9 GeV. 



high AA contribution is an artifact of the ad-hoc off-shell width used. Since the numerical difference 
in varied vector off-shell widths is not that high, the numbers for VV are closer. 
The numbers in Eq.(82) are translated into the following branching ratios 

BRyv- (r ^ 7r7i/^) = 1.05 • 10~^ , BR^a (r ^ 7r7i/^) = 0.15 • 10"^ (83) 
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We can also compare the VV value with the narrow-width estimate. Taking into account the 
lowest-lying resonance /o(770) we get 

BRyy (r — > TT'yi'r) ~ BR{t — ^ pi^r) X BR{p — ^ 7r7) ~ BR{t — > 7r~'K^i^r)BR{p — > 7r7) 

~ 25.52% X 4.5 • 10"'' = 1.15 • 10"^ , (84) 

which is quite a good approximation. 

In Table I we display (for two different values of the photon energy cut-off) how the different 
parts contribute to the total rate. For a low-energy cut-off the most of the rate comes from IB 
while for a higher-energy one the SD parts (particularly the VV contribution) gain importance. 
While the VA contribution is always negligible, the IB — V, IB — A and the remaining SD parts 
{VV and AA) have some relevance for a higher-energy cut-off. 





Xo = 


0.0565 


Xo = 


0.45 


IB 


13.09 


• 10-3 


1.48- 


10-3 


IB^V 


0.02 


10-3 


0.04- 


10-3 


IB^A 


0.34 


10-3 


0.29- 


10-3 


VV 


0.99 


10-3 


0.73- 


10-3 


VA 







0.02- 


10-3 


AA 


0.15 


10-3 


0.14- 


10-3 


ALL 


14.59 


• 10-3 


2.70- 


10-3 



TABLE I: Contributions of the different parts to the total rate in the decay r — tt (in units of 
rr^Ti-iy), using two different cut-offs for the photon energy: — 50 MeV {xq — 0.0565) and = 400 MeV 
(xo - 0.45). 

In Fig. 9 we show the pion-photon invariant spectrum. We find a much better signal of the 
SD contributions as compared with the photon spectrum in the previous Fig. 8, which has already 
been noticed in Ref.[4]. Then, the pion-photon spectrum is better suited to study the SD effects. 
In this case, the VA part is identically zero, since this interference vanishes in the invariant mass 
spectrum after integration over the other kinematic variable. Of course, in the VV spectrum we 
see the shape of the p contribution neatly, as it is shown in Fig. 9, where on the contrary the ai 
exchange in AA has a softer and broader effect. The IB — SD radiation near the ai is dominated 
by IB — A, which gives the positive contribution to the decay rate. While near the energy region 
of the p resonance, we find the IB — SD contribution to be negative as driven by IB — V there. 
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In the whole spectrum only the p resonance manifests as a peak and one can barely see the signal 
of ai, mainly due to its broad width. 

C. Results with the resonance contributions for the kaon channel 

Next we turn to the r~ — )■ K~^Vr channel. In this case, there are several sources of uncertainty 
that make our prediction less controlled than in the — t- ii^^Ur case. We comment them in turn. 

Concerning the vector form factor contribution, there is no uncertainty associated with the 
off-shell widths of the vector resonances , that are implemented as we explained before. However 
we have observed that the VV contribution to the decay rate is much larger (up to one order 
of magnitude, even for a low-energy cut-off) than the IB one for C4 ~ —0.07, a feature that is 
unexpected. In this case, one would also see a prominent bump in the spectrum, contrary to 
the typical monotonous fall driven by the IB term. For smaller values of |c4| this bump reduces 
its magnitude and finally disappears. One should also not forget that the inclusion of a second 
multiplet of resonances may vary this conclusion. 

The uncertainty in the axial-vector form factors is twofold: on one side there is a broad band 
of allowed values for 6 a-, as discussed at the beginning of this section. On the other hand, since we 
have not performed the analyses of the decay r — )• KttttUj- modes yet, we do not have an off-shell 
width derived from a Lagrangian for the Kia resonances. In the r — t- Sttz^t- decays, T^^ has the 
starring role [12]. Since the Kia meson widths are much smaller (90 it 20 MeV and 174 it 13 MeV, 
for the Ki{\21{)) and -fCi(1400), respectively) and they are hardly close to the on-shell condition, 
a rigorous description of the width is not an unavoidable ingredient for a reasonable estimate. We 
decided to use the expression inspired by the p width from Ref. [39] . The explicit forms of the Kia 
off-shell widths which we follow are given in Appendix. A. 

Considering all the sources of uncertainty commented, we will content ourselves with giving 
our predictions for the case of C4 = , — 0.07 and |0a| = 58° , 37°. We first show the results of the 
decay rates in Table. II. The lesson we can learn from the numbers in Table II is that the decay 
rate is sensitive to the value of C4, while different choices of 9a barely influence the final results. In 
order to illustrate our predictions, we present the analogous plots to those as we discussed in the 
r~ — )■ TT^ji/r channel for the case C4 = and = 37° , 58° in Figs. 10 and 11. For C4 = —0.07, we 
give the plots in Figs. 12 and 13, where one can see that the VV contribution from the SD parts 
overwhelmingly dominates almost the whole spectrum. 



31 





xo = 0.0565, C4 = -0.07 
\eA\ = 58°(37°) 


Xo = 0.0565 , C4 
\eA\ = 58°(37°) 


Xo = 0.45, C4 ^ -0.07 
=58°(37°) 


Xo = 0.45 , C4 = 
\eA\ = 58°(37°) 


IB 


3.64-10-3 


3.64- 10-3 


0.31 - 10-3 


0.31 - 10-3 


IB -V 


0.69 - 10-3 


0.10- 10-3 


0.83-10-3 


0.12- 10-3 


IB - A 


0.22(0.25) - 10-3 


0.22(0.25) - 10-3 


0.15(0.18) - 10-3 


0.15(0.18) - 10-3 


VV 


58.55-10-3 


1.30-10-3 


29.04- 10-3 


0.66- 10-3 


VA 


^ 0(^ 0) 


^ 0(^ 0) 


0.09(0.09) - 10-3 


0.01(0.01) - 10-3 


AA 


0.13(0.16) - 10-3 


0.13(0.16) - 10-3 


0.12(0.15) - 10-3 


0.12(0.15) - 10-3 


ALL 


63.23(63.29) - 10-3 


5.39(5.45) - 10-3 


30.54(30.60) - 10-3 


1.37(1.43) - 10-3 



TABLE II: Contributions of the different parts to the total rate in the decay r- K^ji^r (in units of 
^T^Ki^), using two different eut-offs for the photon energy: = 50 MeV {xo = 0.0565) and E-y = 400 MeV 
(xo = 0.45) and also different values of the resonance couphngs. The numbers inside the parentheses denote 
the corresponding results with |^a| = 37°, while the other numbers are obtained with = 58°. 

VIII. CONCLUSIONS 

In this article we have studied the radiative one-meson decays of the r: r" — )• {■K/K)~^Vr. We 
have computed the relevant form factors for both channels and obtained the asymptotic conditions 
on the couplings imposed by the high-energy behaviour of these form factors, dictated by QCD. 
The relations that we have found here are compatible with those found in other phenomenological 
applications of the theory. 

One of our motivations to examine these processes is that they have not been detected yet, 
according to naive estimates or to Breit-Wigner parametrizations. We have checked the existing 
computations for the IB part. Adding to it the WZW contribution, that is the LO contribution 
in coming from the QCD anomaly, we have estimated the model independent contribution 

to both decays, that could be taken as a lower bound. The values that we obtain for the vr channel 
are at least one order of magnitude above the already-observed 3K decay channel even for a high- 
energy cut-off on the photon energy. In the K channel, the model independent contribution gives 
a BR larger than that of the 3K decay channel, as well. Only imposing a large cut-off on 
one could understand that the latter mode has not been detected so far. We expect, then, that 
upcoming measurements at B and tau-charm factories will bring the discovery of these tau decay 
modes in the near future. 
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FIG. 10: Differential decay width ( in units of Tr^Kiyr ) of ^lie process K^^v^. including all 

contributions as a function of x with C4 = 0. The top right plot is the close-up of the top left one in region 
of X > 0.5. The bottom left and right plots are to display the compositions of the interference and structure 
dependent parts, respectively. 



We do not have any free parameter in the r~ — t- tt'^Uj- decay and that allows us to make a 
complete study. Since the IB contribution dominates, it will require some statistics to study the 
SD effects. In this sense, the analysis of the vr-photon spectrum (t-spectrum) is more promising 
than that of the pure photon spectrum (x-spectrum), as we have shown. We are eager to see 
whether the discovery of this mode confirms our findings, since we believe that the uncertainties 
of our study are small for this channel. 

As expected, the higher mass of the Kaon makes easier the observation of SD effects. However, 
there are several sources of uncertainty in the t~ — )• K~^Vr decay that prevent us from having 
a definitive prediction for this channel. The most important one either rises some doubts about 
the value of C4, a parameter describing the SU{3) breaking effect, obtained in Ref. [11] or on the 
sufficiency of one multiplet of vector resonances to describe this decay. We point out that the 
inclusion of the heavier multiplets of resonances will not only directly give contributions to the 
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spectra and the decay rates, but also influence the final results in an indirect way by entering 
the resonance parameters given in Eq.(80). As we have shown, the value of this coupling affects 
drastically the strength of the VV (and thus the whole SD) contribution. Besides, there is an 
uncertainty associated with the broad band of allowed values for 9 a. However the A A contribution 
is anyway not important with respect to that on C4. Even smaller is the error associated with 
the off-shell width behaviour of the axial- vector neutral resonance with strangeness, Kih^l- Since 
we have not calculated the relevant three-meson decay of the tau, we do not have this expression 
within RxT yet. We took a simple parametrization including the on-shell cuts corresponding to 
the decay chains Kiu i — {pK/K*Ti). Since the effect of C4 is so large, we expect that once this 
decay mode is discovered we will be able to bound this coupling. 

As an application of this paper, we are working out [7S] the consequences of our 
study in lepton universality tests through the ratios F (t~ —)• 71^1/^7) /F (vr" —)• ^~f^7) and 
r (r" —7- K~Vt^) /r [K^ — )• fi^i'^j) that were considered in Refs. [79-81] in different frameworks. 
The ratio between the decays in the denominators within was studied in [82, 83] and the 

radiative pion decay has been investigated within RxT in [;>()] recently. 
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The off-shell widths of the resonances used in our numerical discussion are taken from Refs. 
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Appendix A: OfT-shell widths of the intermediate resonances 



[10, 12, 39]. The off-shell width of the p(770) is 




sMp 
967rF2 



(Al) 
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where 

(^pq{s) = ^ {mp + niQ)'^] [s - {mp - rnqf] 6 [s - {mp + rnqY 

9{x) is the standard unit step function. 
The off-shell width of the i<'*(892) is 



Tk'{s) 



sMk* 
1287rF2 



while that of the Kia is 



The ai(1260) off-shell width is 



(A2) 



(A3) 



(A4) 



(A5) 



where 



ai 



-s 



192(27r)3F2A4, \Q 



Ml 



(A6) 



Here r^^(Q2) recalls the three pion contributions and r^'^''^(Q^) collects the contributions of the 
KKtt channels. In Eq. (A6) the symmetry factor S = l/n\ recalls the case with n identical particles 
in the final state. The explicit expressions for Tj'^'^^-A* can be found in Ref.[12]. We stress that 
the on-shell width, T^^{M?), is a prediction and not a free parameter. 



Appendix B: Numerical inputs 

In the numerical discussion, unless a specific statement is given, we use the values given in 
Ref.[3N]. For the pion, kaon and K*{892), we use the masses of the charged particles throughout, 
i.e., 

m^ = 139.6MeV, mx = 493.7MeV, Mk* = 891.7MeV . (Bl) 

For the mass of ai(1260), we use the result of M^^ = 1120 MeV from Ref.[12], which has taken the 
off-shell width effect into account. 

For the resonance couplings, once we have the value of the pion decay constant F in the chiral 
limit, we can determine all of the others except C4 through Eq.(80). For the value of F, we use 
F = 90 MeV. The physical pion and kaon decay constants we are using are 

F^ = 92.4 MeV , Fx = 113.0 MeV . (B2) 
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FIG. 11: Differential decay width ( in units of Fr^/fi/^ ) of the process K^^Vr including all 

contributions as a function of \/t with C4 = 0. The top right plot is the close-up of the top left one in 
region of \/t > 1.1 GcV. The middle left and right plots are to display the compositions of the interference 
and structure dependent parts, respectively. The bottom plot is the close up of the middle right one in the 
region of \/i > 1.1 GeV. 



39 






INT(e,=58") 




INT(e^=37") 




IB-V 




IB-A(e^=58") 




IB-A(e^=37") 





FIG. 12: Differential decay width ( in units of Tt^ku^ ) of the process r~ — ;> K~^Vr including all 
contributions as a function of x with C4 = —0.07. The top right plot is to show the compositions of the 
interference parts. The bottom left one is to show the V A and AA contributions. 
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FIG. 13: Differential decay width ( in units of Tr-yKiy^ ) of tlie process K^^Vr including all 

contributions as a function of ^/t with C4 ~ —0.07. The top right plot is the close-up of the left one in the 
region of ^/t > 1.1 GeV. The bottom left one is to show the compositions of the interference parts. 



